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1.  intrftdMgilpn- 


"Locking”  Is  a  phenonenon  associated  with  the  ntinerical  approximation  of 
certain  problems  whose  mathematical  formulations  involve  a  parameter 
dependency.  The  problem  we  consider  here  is  the  analysis  of  el^lstic 
materials  with  the  parameter  being  v.  the  Poisson  ratio.  For  v  close  to 
0.5  (i.e. ,  when  the  material  is  nearly  Incompressible),  it  is  well  known  that 
various  finite  element  schemes  (for  example,  the  h-version  using  piecewise 
linear  polynomials  on  a  triangular  mesh)  result  in  poor  observed  convergence 
rates  in  the  displacements,  for  practical  ranges  of  the  discretization.  This 
is  due  to  locking  (called  dilatation  or  Poisson  locking  in  engineering).  It 
occurs  because  for  the  limiting  caise  v  =  0.5,  the  exact  solution  u  must 
satisfy  the  constraint 

(1.1)  divu  »  0. 

The  imposition  of  (1.1)  on  the  approximation  as  well  is  what  leads  to  locking 
in  this  example. 

There  are  several  other  problems  where  similar  locking  effects  may  be 
observed  —  for  example,  in  plate  and  shell  models,  where  "shear"  and 
"membrane"  locking  occur  when  the  thickness  "t"  is  very  small  amd  in  heat 
transfer  through  anisotropic  materials  where  locking  occurs  when  the  ratio  of 
conductivities  in  different  directions  is  close  to  zero.  For  problems 
involving  locking,  see  [4]  and  Section  53  of  [15]. 

Various  methods  have  beer  suggested  overcome  the  effects  of  locking. 
One  possibility  is  the  use  of  mixed  methods,  which  involve  reformulating  the 
problem  in  a  special  way.  Examples  of  mixed  methods  that  have  been  suggested 
to  overcome  Poisson  ratio  locking  maiy  be  found  in  [11],  aunong  others.  An 
advantage  of  these  methods  is  that  they  generally  yield  good  approximations 
to  the  "pressures"  as  well. 
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Here,  we  shall  concentrate  solely  on  Poisson  locking  i^en  the  accuracy 
of  the  dlsplaceaents  and  energy  (and  not  the  pressures)  is  of  interest.  Our 
goal  is  to  investigate  the  robustness  of  several  finite  element 
approximations  using  the  standard  (as  opposed  to  nixed)  formulation  (also 
called  the  displacement  formulation).  By  a  robust  scheme,  we  mean  one  which 
leads  to  acceptable  error  levels  using  a  practical  range  of  discretization, 
no  nmitter  how  close  the  parameter  is  to  its  limiting  value.  The  use  of  the 
standard  formulation  avoids  the  special  reformulations  required  by  mixed 
methods  and  in  practice  could  be  the  only  one  available  in  the  context  of 
various  commercial  codes.  Hence  it  is  particularly  useful  to  investigate 
the  Eissociated  locking  and  robustness  properties.  In  this  coxmection,  the 
accurate  recovery  of  the  pressures  may  be  accomplished  through  various 
post -processing  techniques  (see  [22],  for  example). 

In  [4],  we  have  developed  a  general  mathematical  theory  for  locking  and 
robustness  and  their  quantitative  assessment.  We  use  this  theory  here  to 
analyze  Poisson  locking.  Accordingly,  in  Section  2,  v^e  prove  some  required 
regulairity  results  and  in  Section  3,  we  adapt  veurious  definitions  suid 
theorems  from  [4]  to  the  problem  at  hand.  In  this  paper,  we  restrict 
ourselves  to  the  case  of  triangular  eind  parallelogram  quasiuniform  meshes. 

The  case  of  curved  elements,  which  is  particularly  Important  in  the  context 
of  the  p-version,  is  discussed  in  [5]. 

Section  5  contains  vsu'ious  locking  and  robustness  results  for  the 
h-version.  In  [19]  (see  also  [20]),  it  was  shown  that  no  locking  results 
when  polynomials  of  degree  p  2  4  are  used  on  triangular  meshes.  We  present 
am  alternate  proof  here  and  also  give /some  results  for  p  <  4.  The  results 

t 

of  [19]  were  restricted  to  triangular  elements.  Here  we  investigate  the  (ise 
of  two  types  of  rectaurgular  elements  as  well  and  show  that  locking  cannot  be 
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avoided  In  either  case,  for  any  p.  Next,  in  Section  6,  we  indicate  how  one 
can  show  optioal  rates  of  convergence  in  the  displacements  for  the  p  and  h-p 
version  uniformly  in  v.  Our  approach  (and  definition  of  locking)  is 
different  from  that  of  [23]  and  [19],  by  which  one  only  gets  optimality  up  to 
an  arbitrary  e  >  0.  Section  7  contains  extensions  of  our  theory  to  general 
3-d  analogs  of  Poisson  locking. 


2.  Regularity  Results. 

2 

Let  QcR,  OeO  bea  bounded,  simply  connected,  polygonal  domain 
M 

with  boiandary  F  ®  ^  F^,  where  F^  are  open  straight  line  segments  with 

i=l 

internal  angles  >  0.  For  S  c  R^,  we  will  denote  by  H*'(S)  the  usual 

Sobolev  spaces  (r  real)  with  II  *11^,  g  and  I'lj.  s  the 

corresponding  norm  and  seminorm  respectively.  For  any  space  V,  ^  will 

denote  V  x  V  (the  norm  of  V  and  ^  will  be  denoted  by  the  same  symbol), 
(k) 

We  will  use  C  (S),  k  i  0  integer  to  denote  the  usual  set  of  functions  o 


S  with  k  continuous  derivatives,  and 


superscript  k  being  omitted  when  0). 


c'^^S) 


to  denote  its  norm  (the 


The  problem  we  are  Interested  in  is  the  elasticity  problem  given  by  the 
following  Lam^-Navier  equations 


(2.1) 


_*•-»  -E  E 

\  *  2n^^%  ■  2T1+pT( 


graddivUy  =  ?  in  fl 


(2.2) 


T^(u)  =  g  on  F 


where  u^  =  (Uj.u^)  and  where  the  tractions  (T^(Uy))j  are  given  for  i 


1.2  by 


(2.3) 


[ife]  I  ('ij‘%’  * 
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Here,  (nj.ng)  Is  the  unit  outuard  nornal  to  r  and  Is  the  strain 

tensor  given  by 


ij 


1  fauj^  8^") 

2  (^9xi  8xjJ  • 


The  coefficient  0  ^  v  <  0.5  represents  the  Poisson  ratio  and  E  the 
modulus  of  elasticity,  which  are  related  to  the  Lam6  constants  by 


(2.4) 


Ev 

(l+w)(l-2v) 


^  “  2(l+v)  • 


We  assume  that 

(2.5)  JJ?-^dx  +  Jg-^ds  »  0 

Q  r 


for  any  rigid  body  motion  ^  to  ensure  that  (2.1)  and  (2.2)  have  a  solution 
(unique  up  to  S). 

As  usual,  the  components  1  1,2  of  the  stress  tensor  are  then 

given  by 


X8*'  + 


o  aul 


\e 


’Suj  ^  8^1 
8xi  8X2J 


1/  ^ 

where  0  »  dlv  u  , 
v 

We  have  used  the  index  v  to  emphasize  the  dependence  on  v.  If  no 
ambiguity  can  occur,  we  will  omit  this  Index. 

We  will  assume,  for  simplicity  (but  without  loss  of  generality),  that  E 
»  1.  We  will  also  consider  the  following  variational  form  instead  of 
(2.1) -(2.2):  Find,  for  given  u  e  [0,  0.5),  a  u^  €  3^(0)  satisfying, 

Vv  c  S^(n), 
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(2.6) 


V 


where 


B^(u^.v)  -  a^(u^.v)  +  dlvv)  =  F{v) 


(2.7) 


(2.8) 


aj^(u,v) 


F(v) 


*  ff  Z 

h 

=  j*j*?*vdx  +  J*g*v  ds. 


We  denote  the  problem  (2.6)  by  P^.  It  can  be  skssumed  equivalent  to 
(2.1) -(2.2)  without  loss  of  generality. 

Corresponding  to  (2.6),  we  define  the  energy  norm  by 

l|3|||  y  =  By(3.u).  ueS^n). 

Using  Korn’s  Inequality,  we  see  that 

(2.9)  C,|3|j  „  S  ilSlj,  ^  S  C2(1-2»)-’^  ISIj  „ 

where  ||‘|l|^  q  denotes  the  norm  In  the  quotient  spaice 

i{^(Q)\{Rlgld  Body  motions)  and  and  C2  are  constants  Independent  of  1/. 
Obviously,  for  y  bounded  away  from  0.5,  the  two  norms  are  equivalent. 

For  euiy  k  2  1,  0  S  y  <  0.5,  let  us  define  the  space  \  y  “  ^ 

furnished  with  the  weighted  norm  (modulo  rigid  body  motions). 


(2.10) 


Related  to  this  definition  aire  the  characterizations  of  balls 


(B  e  R  )  given  by 


»  (3  €  a^(n),  S  B) 


The  speu:es  ^  are  the  natural  spaces  to  consider  while 
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characterizing  the  solutions  of  (2.6)  (or  (2. 1)  -  (2.2)).  Their  choice  is 
motivated  by  the  following  theorea  from  [23]. 


Theorem  2.1.  Let  Q  be  a  smooth  domain  with  smooth  boundary  T.  Let  k  be 
an  integer  21.  If  u^  6  3^(0)  denotes  the  solution  to  (2.1) -(2.2)  for 
data  ?  6  ^(Q),  g  €  ^^(F),  then  u^  €  1^  ^  with 

B  =  Q  ♦  "8!tjj_3/2,r) 

where  C  depends  on  k  but  is  independent  of  v. 

For  the  case  that  the  domain  is  a  polygon,  the  above  theorem  will  again 
be  valid,  but  only  for  a  restricted  range  of  k,  i.e.,  k  S  kjj,  where  k^ 
is  determined  by  the  domain.  We  now  prove 

Lemma  2. 1.  Let  u^  €  ^  for  k  2  2,  0  S  u  <  0.5.  Let  ?  and  g  be 
defined  in  terms  of  u^  by  (2.1),  (2.2)  respectively.  Then  ?  and  g 
satisfy  (2.5)  and 

'?'k-2,n  *  Z  >*'ic-3/2.r,  ^  “ 

i 


where  C  is  a  constant  Independent  of  v,  u^  and  B. 

Proof.  Letting  v  in  (2.6)  be  a  rigid  body  motion,  it  is  easy  to  see  that 
(2.5)  is  satisfied.  From  (2.1),  we  see  that  (E  »  1) 

’*^®k-2,Q  ^  2(l+v)  **^%**k-2,n  *  2il*v)[l-2v)  “i;”k-2,n 


^  2  [*’%**k,n  1-21/ ’’‘*^''%**k-l,n] 


Similarly,  by  (2.2),  (2.3), 
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^  **®k-3/2,ri  ^  1+y  ^  [’*“y*'k-l/2,ri  *  l-2w  ^••k-3/2.  rj 


S  Cgu  II.  S  CB 


which  gives  (2.11). 

We  will  be  interested  in  the  limiting  sets 

=  (u  e  y  Vy  e  10,  0.5)}  . 

Obviously,  these  may  be  equivalently  characterized  as 


where 


=  (3  €  rf'(Q).  ||3||j^  S  B> 


Hj^  =  {u  €  ii^'cn),  dlvu  = 


0} 


Let  us  prove  the  following  Lemma. 


Lemma  2.^.  Let  u^  €  1^  k  k  2,  0  S  y  <  0.5  be  such  that  is  satisfies 
(2.1)  with  ?  *  3.  Then  there  exists  a  u^^  €  such  that 

-  \'k.n  ^  ’'^"-2'” 


(2.12) 


where  the  constaoits  a  and  y  are  independent  of  y,  u^^  and  B. 

Proof.  Since  u^  satisfies  (2.1)  with  ?  =  0,  the  stresses  cr^j  relate 
to  the  Airy  biharmonic  function  U  in  the  usual  fashion, 

(2.13) 


a^u 


d^\j 


a^u 


11  a  2’ 
3X2 


12  3xi3x2’ 


22  .2 
3X1 


Let  us  denote  (up  to  an  arbitrary  constant) 

(2.14)  P  =  AU. 

Then  since  U  is  biharmonic,  P  is  harmonic.  Denoting  by  Q  the  harmonic 
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conjugate  of  P  and  using  the  fact  that  Q  is  simply  connected,  we  have  for 
any  s  ^  1, 


(2. 15) 


Also,  with  2  =  Xj  +  iXg, 


''“'s.a- 


h(z)  =  PCXj.Xg)  +  iQCXj.Xg) 


(2.16) 

is  a  holofflorphic  function  on  Q.  We  define 


(2. 17)  ^(z)  **  J  dz  =  *  ^P2^^r^2^  ‘ 

Then  (see  eq.  (30.8)  of  (171),  the  solution  to  (2.1) -(2.2)  is  given  (up  to 
rigid  body  motion)  by 

(2. 18) 


where 

(2.19) 


“i  ■ 


Cj(y)  »  -d+u),  C2(y)  *  4(l+p)(l-y) 


By  equation  (30.7)  of  [17],  we  have 

dxi  4 

so  that  using  (2.18)  and  (2.14), 

dlv 


1  =  1,2 


+ 


%  -  (s'">  *  "^1'’ 


*  (l-2y)(l+y)P  . 


Hence, 


'’^’*k-l,n  ^  (l-2u)(l+y)  ’*'*^''%"k-l,n  ^  "%"k,y  ^ 

Therefore,  using  (2. 15)  -  (2. 16) ,  for  k  >  2, 

llhllk.i^jj  "  CB 

(where  the  norm  is  a  quotient  norm  modulo  constants). 
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Then,  by  (2.17),  for  i  •  1,2  (modulo  rigid  body  motions). 


(2.20) 

so  that  by  (2. 18)  -  (2. 20) , 

(2.21) 


'Pl»k,Q  ^ 


5  ^  •iPi^k.n 


^  CB  . 

We  now  define  (up  to  rigid  body  motion)  by 

(2.22)  u[  =  Cj(O.S)  +  C2(0.5)  pj,  i  =  1.2. 

Then  dlvu^^  =  0  and  by  (2.20)  -  (2.22), 

"\”k,Q  " 

so  that  u^  *  (2.18),  (2.22), 

(2.23)  uj-u^.  (l-2v)[|||.*  (U2p)p,] 

so  that  by  (2.20),  (2.21),  (2.23), 

'“a  -  \'k.O  ^  »B(l-2y), 

which  establishes  the  lemma.  □ 

We  now  extend  the  above  lemma  to  the  case  that  ?  ^  5,  to  get  the 
following  theorem. 

Theorem  2.2.  Given  u^  e  k  2:  2,  0  S  r  <  0.5,  there  exists  a  Uj^  e 

/  ayp\ 

^  such  that  (2.12)  holds  with  a, 7  being  constamts  independent  of  v, 
and  B. 
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t3 


Proof.  Let  e  k  a  2  be  given.  Define  ?  by  (2.1).  We  will 

reduce  this  case  to  the  case  of  Lemma  2.2,  where  ?  -  5. 

Vte  first  find  a  particular  solution  of  (2.1).  Since  by  Lemma  2.1,  ?  e 

i?'  ^(Q)  for  some  k  2  2,  we  may  find  an  extension  ?  of  ?  to  the  whole 
of  ^  such  that  ?  has  compact  support,  satisfies 

(2.24)  4  C||?||^.3_„  S  CXB 

(by  Lemma  2.1)  and  also  satisfies  the  compatibility  condition 

(2.25)  JJ  ?  •  S  a  0 

for  any  rigid  body  motion  ft.  Then  the  problem  (2.1)  with  ?  replaced  by  ? 

2 

will  have  a  unique  solution  over  R  (up  to  rigid  body  motion),  which  we 
denote  by  Obviously,  is  a  particular  solution  for  (2.1). 

Now,  let  W„  , .  i  «  1,2  denote  the  Fourier  transform  of  (ft  ),  and 

1/  f  1  V  1 

the  Fourier  treuisform  of  F^.  Then  if  represent  the  trEuisformed 

variables,  we  have  by  (2.1), 


*  1-24* 

4Vl*«.<2V2]]-f 

2nk((^*4]v2*  1^4* 

iK.2  *  «l«2“».l])  - 

which  gives 
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Using  (2.26),  we  see  that 


(2.28) 


div  «A.2^ 


(l-i;)(C‘+€‘) 


and  putting  t/  =  0.5  in  (2.28),  we  obtain 


divWj^  *  0. 


(2.29) 

Now  by  Parseval’s  equality, 

(2.30)  !»}  I,  _  i  IvJ  I,  _2  »  |i^  I,  _2  . 

wl,n  u*l,R‘‘  ^'1,8** 

Using  (2.26),  we  have 

2 

'  =  ^  I  '«* 

1«1 


s  c 


C|?I^.2,r2 


1=1 


5  C(tB) 

by  (2.24).  Hence,  by  (2.30)  -  (2.31), 


(2.32) 

Similarly, 

(2.33) 


||W  II.  _  i  CtB  . 

"  p"k,n 


«\Nk,n  ^ 


Also,  by  considering  (2.28),  we  obtain 


(2.34) 


PdivWy||j^_j  6  Ct  (l-2y)B 
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so  that  conblnlng  (2.32)  and  (2.34)  gives 


(2.35)  iJ  € 

u  K,y 

for  some  <c.  Moreover,  by  (2.29)  and  (2.33), 

(2.36) 

Next,  using  (2.26)  smd  (2.27),  vie  have 


w  ,-w,  , 
v,l  L,  1 

l-2w 

A 

h 

W 

v,2  L,2 

(l«)(?”+^)-25= 

.*^2 

so  that  using  an  argument  similar  to  (2.31),  we  have  for  some  t, 

(2.37)  ^  fB(l-2u)  . 


Now  let  on 


(2.38) 


Then  by  (2.35),  w^  e  foi*  some  a.  Moreover,  w^  satisfies  (2.1) 

?  =  5.  Hence,  applying  Lemma  2.2,  there  exists  a  w^^  €  such  that 


(2.39) 


II Wy  -  Q  ^  rB(l-2i;). 


Finally,  taking 


(2.40) 


we  see  (using  (2.36)) 


„(aB) 

"k.L 


Also  by  (2.38),  (2.40) 


"<«„  *  -  ‘«L  *  “L>“k.n 


s  ||W  -  W,  L  „ 
V  L"k,n 


w  - 

V 


L"k,n 


i  rBil-2v) 


with 
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for  some  7,  using  (2.37),  (2.39).  This  proves  the  theorem. 


□ 


3.  Locking  and  robustness. 

We  now  discuss  the  approximate  solution  of  (2.6).  Let  {^}  be  a 
sequence  of  finite  dimensional  subspaces  of  il^(Q)f  (N  denoting  the 
dimension,  N  €  /().  Vfe  then  find  u|^  €  ^  satisfying 

(3.1)  By(Uy,v)  ■  B^(Uy,v)  Vv  €  ^. 


(3.1)  immediately  gives 

(3.2)  llUy 


E.v 


inf  II U  -  w|| 


E.V 


The  sequence  {^}  defines  an  extension  procedure  9,  l.e. ,  a  rule  to 
Increase  the  dimension  N  (and  thereby  decrease  the  error  in  (3.2)). 

We  will  restrict  our  attention  to  the  case  when  the  exact  solutions 
belong  to  the  sets  H.  c  3^(0)  (k  2  2)  Introduced  in  the  previous 
section.  We  assume  that  the  sequence  9  -  is  such  that  for  ainy 

0  S  V  <  0.  S, 


(3.3) 


CjFq(N)  S  Sup 


-  ^'i.n  *  Vo"" 


where  Fq(N)— >0  as  N— ♦«,  Fq  independent  of  v,  and  C^,  Cg  independent 
of  N  and  v. 

Let  0  <  Vq  <  0.5  be  bounded  away  from  0.5.  Then  using  (2.9),  (3.2) 
and  (3.3),  we  see  that  the  following  will  hold  uniformly  for  all  0  ^  y  < 

(3.4)  Cj(yQ)FQ(N)  5  Sup  £^(3^  -  3”)  S  C2(yQ)FQ(N) 

u  eH? 
y  Tc,y 

where 

(3.5)  E^(w)  *  llwllj  Q  or 
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and  C^,  now  depend  on 

A  procedure  9  for  vdiich  the  estimate  (3.4)  holds  xinlformly  for  all 
0  ^  r  <  0.5  will  be  called  free  from  locking  for  the  sets  ^  with 
respect  to  the  measure.  We  make  this  more  precise  by  using  definition 

3.1  below,  which  has  been  adapted  from  [4],  in  which  a  more  general  treatment 
may  be  found  (for  e.g. .  we  could  formulate  the  question  of  locking  in  terms 
of  other  error  measures  and  solution  sets  different  from  the  ones  considered 
here ) . 

For  V  €  [0,  0.5)  and  N  €  we  define  the  locking  ratio  with  respect 
to  the  spaces  ^  c  i^'dl)  and  error  measures  {E^>  (as  in  (3.5))  for  the 
problems  (3.1)  by 


(3.6) 


L(i;,N) 


Sup  E„(S_  -  S^XFglN)) 

U  €H? 

V  K,V 


We  then  have 


Definition  3. 1.  The  extension  procedure  9  is  free  from  locking  for  the 
family  of  problems  (3.1),  v  e  [0,  0.5)  with  respect  to  the  solution  sets 
y  c  3^(0)  euid  error  measures  E^  if  and  only  if 

Llm  sup  sup  L(r, N)  =  M  <  « . 

N— »  »  Ve(0,  0.5)  J 

9  shows  locking  of  order  f(N)  if  and  only  if 


where  f(N)— as 
most)  order  f(N) 


0  <  Lim  sup  jsup  L(r,N) 

N — >  00  ^  w  ' 


C  <  oo 


N— )oo.  It  shows  locking  of  at  least  (respectively  at 
if  C  >  0  (respectively  C  <  oo). 


Related  to  the  concept  of  locking  is  the  concept  of  robustness,  which  we 
define  as  follows: 
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Definition  3.2.  The  extension  procedure  f  Is  robust  for  the  feunlly  of 
problems  (3.1),  v  €  [0,  0.5)  with  respect  to  the  solution  sets  it  c 

K,  V 

a^cn)  and  error  measures  If  and  only  If 

Llm  Sup  Sup  E  (u  -  u.  )  ■  0  . 
u  u  €hJ 

V  ic.v 

It  Is  robust  with  \inlform  order  g(N)  If  and  only  If 

Sup  Sup  ^  Ey(u^  -  uJJ)  S  g(N) 
p  u  sH, 

W  K,U 

where  g(N)— ^0  as  N^oo. 

The  relationship  between  the  above  two  definitions  Is  given  by  the 
following  theorem,  from  [4]. 

Theorem  3.1.  f  is  free  from  locking  if  and  only  if  It  Is  robust  with 
uniform  order  F^CN).  Moreover,  let  f(N)  be  such  that 

f(N)FQ(N)  =  g(N)-»0  as  N-)». 

Then  f  shows  locking  of  order  f(N)  if  Buid  only  If  it  is  robust  with 
uniform  order  g(N). 

Note  that  9  Is  non-robust  If  and  only  If  It  shows  locking  of  order 

(Fq(N))"^ 

Let  us  make  some  comments  about  the  meaning  of  the  notions  we  have 
Introduced.  The  Inequality  (3.3)  chareicterlzes  the  approximation  properties 
of  the  space  ^  with  respect  to  the  set  i.e.,  the  smallest  error 

which  could  be  achieved  (for  the  most  unfavorable  w  in  i^).  The  locking 
ratio  characterizes  the  ratio  of  the  accuracy  of  the  finite  element  solution 
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(for  the  most  unfavorable  exact  solution  in  H.  )  to  the  best  accuracy 

K  f 

which  could  be  achieved  In  the  sense  of  (3.3),  using  the  finite  element  space 
V^.  Hence  an  extension  procedure  (used  to  construct  a  sequence  of  finite 
element  solutions)  is  called  free  of  locking  if  for  all  u^  €  v  e.  [0, 

0.5),  it  leads  to  a  rate  of  convergence  which  is  asymptotically  the  same  as 
the  best  rate  vdiich  the  set  could  provide  (in  the  sense  of  (3.3)).  The 
order  f(N)  of  locking  characterizes  the  asymptotic  strength  of  the  locking, 
i.e.  the  convergence  rate  of  the  finite  element  solution  compared  to  the  best 
possible  rate.  The  robustness  gives  a  measure  of  the  rate  of  convergence 
which  holds  independent  of  the  parameter  v. 


1  ■  In  [19],  [20],  [23],  the  condition  used  to  obtain  uniform  rates 
of  convergence  in  v  (i.e.,  characterize  the  absence  of  locking)  is 


(3.7) 


1%  -  “U'l.D  ^  -  ’<1.0  ‘  ">•  O- 


where  C  is  a  constant  independent  of  N,  v.  Here,  u^,  u”  are  solutions 
corresponding  to  fixed  ?  and  g,  with  different  v.  (3.7)  is  a  stronger 
condition  than  the  one  in  Definition  3.1  —  it  is  equivalent  to  the  so-called 
"divergence  stability"  and  is  necessary  EUid  sufficient  for  the  absence  of 
locking  when  both  the  displacements  and  pressures  EU'e  considered  together. 

It  can  be  characterized  using  Definition  3. 1  by  choosing 

e„(3)  -  '‘“’""o.o 


in  the  locking  ratio  (3.6).  Alternately,  one  could  re-define  the  locking 


ratio  (3.6)  as 


L(i;.N) 


E  (u  -  u  ) 

V  V  V 

-  . 

inf  E  (u  -v) 

V  V 
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and  then  use  Definition  3.1  (with  E^(w)  »  Uw||^  q). 

As  mentioned  in  the  introduction,  we  restrict  our  investigation  here  to 
the  displaiceaents  alone,  with  Definition  3.1  (and  not  (3.7))  chairsicterizing 
the  absence  of  locking.  This  leads  to  a  relaxation  of  a  condition  needed  In 
a  theorem  from  [19]  (see  Theorem  5.4)  and  an  improvement  in  the  robustness 
result  from  [23]  (see  Section  6).  □ 

As  analyzed  in  [4],  the  question  of  locking  meiy  be  reduced  to  one  of 


approximability  alone  if  the  problem  and  solution  sets  satisfy  a  certain 
condition,  called  condition  (a).  This  condition  requires  that  for  any  u^  e 
y,  there  exist  a  Uq  c  (for  some  B'  independent  of  u^,  r;  u^ 

depending  on  u^)  such  that 


'“u  -  “o"k.O  * 


with  C  a  constant  independent  of  i;  and  u^.  Essentially,  condition  (a) 
says  that  the  solutions  u^  are  close  enough  to  functions  u^  in  the  limit 
space.  As  shown  in  [4],  this  allows  us  to  answer  the  question  of  locking  for 
solutions  u^  Just  by  considering  the  approximability  of  functions  in  the 
limit  space  and  lesids  to  some  useful  theorems  concerning  locking  euid 
robustness. 


For  our  particular  problem,  we  see  that  condition  (a)  is  satisfied  (with 
B'  =  oB)  by  equation  (2.12)  (Theorem  2.2).  Then  Theorem  2.4  from  [4]  holds 
and  can  be  stated  eis  follows. 


Theorem  3.2.  Let  us  consider  the  family  of  problems  (3.1),  v  €  [0,  0.5) 
with  the  solutions  sets  „  C  ifcn).  k  k  2.  Then  the  extension  procedure 
y  is  free  from  locking  with  respect  to  the  S^(n)  norm  if  and  only  if  it  is 
free  with  respect  to  the  energy  norm.  It  shows  locking  of  order  f(N)  in 
the  3^(0)  norm  if  and  only  if  it  shows  locking  of  order  f(N)  in  the  energy 
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norn. 


Theorem  3.2  shows  that  the  locking  behavior  using  the  two  error  measures 
under  consideration  is  the  same.  This  is  important,  because  for  theoretical 
purposes,  the  2^(0)  norm  is  easier  to  work  with,  while  computationally, 
results  obtained  by  programs  like  MSC/PROBE  are  often  in  the  energy  norm  (see 
[5]).  We  shall  therefore  only  refer  to  the  locking  of  ?,  without 
specifying  an  error  measure. 

Let  us  define,  for  any  space 

Z(^)  =  {u  €  ^,  dlvu  =  0>. 

Then  we  may  reduce  the  question  of  locking  to  one  involving  approximability 
alone. 

Theorem  3.3.  For  the  problems  (3.1)  with  solution  sets  H.  ,  k  i  2,  the 
extension  procedure  f  is  robust  with  uniform  order  g(N)  given  by 

(3.8)  g(N)  «  Sup  inf  1|2  -  w||  « 

^  3€Z(^") 

Also,  with  (^q(N}  as  in  (3.3),  9  is  free  from  locking  if  and  only  if 

(3.9)  g(N)  S  CFq(M). 

It  shows  locking  of  order  f(N)  if  aind  only  if 

CjFQ(N)f(N)  5  g(N)  ^  C2FQ(N)f(N). 

Proof.  The  theorem  follows  using  (2.12)  euid  Theorem  2.2(B)  of  [4].  □ 

We  see  therefore  that  in  order  to  check  for  locking,  we  need  only 
estimate  g(N).  We  now  derive  an  equivalent,  more  convenient  formulation  of 
g(N). 

Let  for  2uiy  space 

W(^)  =  f  Curl  #  e  ^>. 
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Then  it  is  easy  to  see  that 


(3.9) 


Z(^)  =  (u  =  Curl^,  ^  c  W(^)}  . 


Lemma  3. 1.  For  kil,  and  only  if  there  exists  ^ 


H  Q  ^  such  that 


(3. 10) 


u  *  Curl ^ 


tdiere  a  is  a  constant,  independent  of  u,  B. 

/^yp\  ^  ^ 

Proof.  Let  ^  e  Hjj+j  (0)-  Then  defining  u  by  (3.10),  we  see  that  u  e 
“k.L 


Hence  u  € 

Conversely,  let  u  €  Then  we  may  find  ^  satisfying  (3.10). 

Since  ^  is  arbitrary  up  to  a  constant,  we  may  choose  it  such  that  [  ^  »  0 

•'q 


so  that 


Then 


k+i,n 


so  that  ^ 


Let  g(N)  be  as  defined  in  (3.8).  Let 


(3. 11) 


Then 


g(N)  *  Sup 


inf  11^  -  x\\. 


z€W(^) 


Cjg(N)  S  g(N)  S  C2g(N) 
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with  C^,  C2  independent  of  N. 

Proof.  Let  u  c  ^  eind  v}  c  Z(^).  Then  by  Lemma  3.1  and  (3.9)  we  can 
find  ^  e  H^J^^(n),  %  €  W(^)  such  that 

(3. 12)  u  *  Curl  w  =  Curl  %  • 

Hence, 

Hu  -  will  Q  *  yCurl^  -  Curljcllj  Q  i  ||^  -  zlig.n 

so  that 


g(N)  S  CggCN). 

Conversely,  let  ^  e  Hj^^^(n),  %  €  W(^).  Then  using  Lemma  3.1,  let  u  6 

we  Z(^)  satisfy  (3.12).  By  adjusting  z  by  a  constamt,  we  can 
ensure  that  (3.12)  remains  true  with 


SO  that 


inf  110  -  zlig  jj 
Z€W(V'‘) 


S  Cinf  IlCurl  (0-z)lL  „ 
X€W(V'') 


*  C  inf  I 

;:eZ(^") 


u  -  w|| 


i,n- 


The  theorem  follows. 


□ 


4.  The  extension  procedure  f. 

We  now  define  the  subspaces  Let  {y^}  be  a  quasiuniform 

sequence  of  meshes  on  (2  (in  the  usual  sense,  see  [13]),  parauneterized  by  h, 
the  mesh  spau:ing.  Each  7^  consists  of  straight -sided  parallelogram  or 
triangular  elements  (1^,  i  s  l,2,...,n(h)  such  that  ^  r\  is  either 
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enpty,  a  vertex  or  an  entire  side  for  i  ^  J-  Vfe  will,  in  particular,  be 
interested  in  the  case  that  R  is  a  domain  which  can  be  covered  by  uniform 
triangular  or  rectangular  meshes  5^,  Jg  of  the  type  shown  (for  a  unit 
square)  in  Figxire  4.1. 


The  general  quasiuniform  versions  of  these  meshes  will  be  denoted  by 
respectively.  The  meaning  of  quasiuniform  triangular  (7^)  meshes  on  a 
polygonal  domain  R  is  clear.  The  quasiuniform  meshes  7^  will  consist  of 
elements  R^  which  are  rectangles  (of  possibly  different  sizes,  and  not  just 
uniform  sqiiares).  We  will  also  consider  briefly  the  case  of  meshes  7^ 
consisting  of  parallelograms. 

For  S  a  triangle  or  parallelogram,  let  Pp(S)  denote  the  set  of 
polynomials  on  S  of  total  degree  i  p.  For  S  a  parallelogram,  we  let 
Q^(S)  be  the  set  of  all  polynomials  on  S  with  degree  ip  in  each 
variable  and  define 

Q|^(S)  =  Pp(S)  ®  Xj)^>, 

to  be  the  space  of  serendipity  elements  (see  [13]). 
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For  any  aesh 


we  now  define 


)} 


eund  for  k  2;  0, 

pII  V  “  =  pI!  -1 

P»  K  p«  K  p»  1 

For  rectangular  meshes  ^  (-4  or  we  analogously  define 

qU  .  (y^)  and  qJ;  '(y^)  with  replaced  by  Q^tQ^),  Q'(n?) 

p,  K  P,  K  pi  plpl 

respectively.  Then  our  finite  element  space  ^  ®  ^  will  be  taken  to  be 

one  of  or  3^  . 

Pi  0  Pi  0  P,  0 

We  now  estimate  in  (3.3). 

Theorem  4. 1.  Let  f  =  {^}  =  (^>  consist  of  the  h-version  using  a  family 
of  qxiasiuniform  meshes  <y^}  of  parallelograms  and  triangles,  with  degree  p 


i  1  fixed.  Then  (3,3)  is  satisfied  with  Fq(N)  given  by 


(4.1) 


Fq(N)  ~  CN 


-p/2 


CN 


-(k-l)/2 


for  p  <  k-1 
for  p  k  k-1 


where  C  is  independent  of  N  but  depends  on  p,k. 


Proof.  We  note  that  by  the  usual  results  (see  (131)  for  the  h-version, 

Fq(N)  S 

-2 

since  N  *  0(h  ).  The  result  for  p  i  k-1  follows  by  noting  that 

Fq(N)  i 

using  the  theory  of  n-wldths  [18].  For  p  <  k-1,  the  corresponding  lower 
bound  has  been  proved  in  Lemma  3.2  of  [4l. 

D 
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Theorem  4.2.  Let  f  *  {^}  =  {^}  consist  of  the  p-version  using  a  fixed 

P 

mesh  with  increasing  polynomial  degrees  p— >od.  Then  (3.3)  holds,  with 

Fq(N)  given  by 

(4.2)  Fq(N)  ~ 

Moreover,  ais  p — »«,  (4.2)  also  holds  if  the  h-p  version  over  a  quaisi- 
vinlform  family  of  meshes  {7^}  is  used. 

Proof.  We  refer  to  [3]  for  the  proof  of  the  p-version  and  [2]  for  the  proof 
of  the  h-p  version.  □ 

Remark  4. 1.  In  the  case  of  the  h-p  version,  the  following  more  refined 
estimate  has  been  established  in  12] 

(4.3)  Fq(N)  - 
where  N  =  N(h,p). 


Let  us  now  chsu'acterlze  the  spaces  W(^)  for  the  various  choices  of 

P 

V^.  These  will  be  used  in  the  next  section. 

P 

First,  we  see  that  over  a  single  element  (T  =  triangle,  S  = 


pEU*allelograra), 

(4.4) 

«Pp<T>]  - 

(4.5) 

“(V®’]  =  Vi'®’ 

U  9p<S> 

(4.6) 

*'(3p<s)]  -  Vi‘=’ 

for  ] 

=  P3(S)  U  Q2(S)  = 
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Therefore, 


(4.7) 


W(?^  )  *  P*' 

Vi.i 


for  the  trlaingular  meshes  5^,  Also,  for  rectaingular  meshes 


(and  yj). 


(4.8) 


“<o> "  ('^*1.1  U  <.] 


and  for  p  *  2 


(4.9) 


W(3^  ')  =  j  . 

p,0  p+1,1 


with  the  case  p  =  2  being  given  by  (4.8), 


5.  The  h-verslon. 

Ue  now  analyze  the  locking  and  robustness  of  some  h-verslon  extension 
procedures.  Theorem  3.4  will  play  a  key  r-^le  in  our  analysis. 

We  first  analyze  triangular  meshes. 


Theorem  5. 1.  Let  the  extension  procedure  9  consist  of  the  h-version  with 
piecewise  polynomials  of  degree  1  using  a  uniform  mesh  3^  as  In  Figxire 
4.1.  Then  with  respect  to  the  solution  set  H.  ,  k  2:  2,  9  shows  locking 

V 

1/2 

of  order  N  ,  i.e.,  the  extension  9  Is  not  robust. 


Proof.  We  have  ^  ■  ^.0  and  by  (4.7)  W(^)  “  (*2  1‘  Theorem  3.4,  we 

have  to  estimate 


g(N)  » 


SUD. 


Assume  that  g(N)— >0  as  N— »aa.  Then  by  the  embedding  theorem  we  also  have 


inf 

xeP 


h 

2,1 


as 


N — >« 


But  by  Theorem  3  of  [7],  the  set  is  dense  in  C(n)  only  if 

s  ^  3.  Hence  we  have  a  contradiction  and  the  theorem  is  proven. 

Let  us  now  define 


a 


r(p) 


p  for  p  =  2,  3,  4 

p  + 1  for  p  S  5. 


Then  we  have  the  following. 


Lemma  5. 1.  Let  be  the  uniform  trianguleu'  mesh  of  Figure  4. 1. 


Then  for 


u  €  H^(n).  p  S  2, 


(5.1) 


inf 

V^ph  (jh) 

p,  1  1 


Moreover,  there  exists  a  function  Q  €  C  (tl)  satisfying 


(5.2) 


inf  IIQ-vIl-  _  ^  Ch' 

veph  (yh)  o,n 

p,  1  1 


r(p) 


Proof.  The  result  for  p  2  5  is  staindard  (see,  for  e.g. ,  [12]).  For 

2  ^  p  £  4,  (5.1)  is  a  generalization  of  Theorem  4  of  [7],  and  has  been 

proven  by  C.  de  Boor  [6].  The  L  analog  of  the  lower  estimate  (5.2)  was 

00 

established  in  [8]  for  p  =  3  and  [9]  for  p  =  4.  These  results  have  been 
recently  generalized  to  the  L^  (1  ^  q  ^  «)  case  in  [10].  □ 

Lemma  5.1  results  in  the  following  lemma.  The  proof  of  (5.4)  is  due  to 
C.  de  Boor,  R.  DeVore  and  A.  Ron. 


Lemma  5.2. 
(5.3) 


For  above,  u  €  8*^(0),  p  ^  2, 

inf  ||u  -  v||  n  S  Ch”^"  (k,  r(p)  )-2||^|| 
veP*'  (?“)  ’ 

p,  1  1 


k,n  • 
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Moreover,  there  exists  a  function  Q  e  c“(n)  satisfying 

(5-4)  inf  llQ-vll 

veP**  (7“) 

p.  1  1 

Proof.  (5.3)  follows  from  (5.1)  by  using  an  inverse  property  argument.  To 
prove  (5.4),  let  us  choose  the  same  function  Q  from  (5.2)  in  Lemma  5.1. 
Suppose  V  €  pIJ  ,  (j!')  is  such  that 

p,  1  1 

Then,  by  using  the  above  and  (5.1),  there  exists  w  €  P^  ,(?^)  such  that 

P.  1  1 

IKQ-v)  -  w||q  Q  i  Ch^llQ-''ll2,Q 

=  o(h^^Ph. 


which  is  a  contradiction  to  (5.2). 

Lemma  5.2  gives  sharp  upper  and  lower  bounds  for 


□ 


g(N)  =  Sup  inf  ll♦-aclL  o 

*€H®  (n)  xeP**  (7“) 

1  P  p 1  1 


for  k  S  r(p)  -  1.  In  fact,  we  obtain 


(5.5) 


C^h*'^P)-2  ^  ^  j,^j^r(p)-2 


Theorems  5.2  -  5.3  below  follow  immediately  from  (5.5). 


Theorem  5.2.  Let  the  extension  procedure  3  consist  of  the  h-verslon  with 

piecewise  polynomials  of  degree  5  2  on  7^.  Then  with  respect  to  the  set 

y,  extension  7  is  uniformly  robust  with  order  N  for  k  >  2  and 

l/p 

shows  locking  of  0(N  )  for  k  i  3. 
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Theorem  5.3.  Let  the  extension  procedure  5  consist  of  the  h-verslon  with 

piecewise  polynomials  of  degree  S3  on  5^.  Then  with  respect  to  the  set 

extension  9  Is  uniformly  robust  with  order  for  k  i  3  and 

1/2 

shows  locking  of  0(M  )  for  k  ^  4. 

In  Theorems  5.1  -  5.3  (for  p  S  3),  we  restricted  ourselves  to  the 
uniform  mesh  7^.  It  was  essential  that  the  mesh  have  three  pairwise 
Independent  directions,  since  with  meshes  with  additional  lndep>endent 
directions,  the  results  will  be  different.  For  p  >  3,  we  may  drop  the 
restriction  of  uniform  meshes  and  prove  the  following  theorem. 

Theorem  5.4.  Let  the  extension  procedure  9  consist  of  the  h-verslon  with 
piecewise  polynomials  of  degree  p  i  4  on  a  general  quaslunlform  mesh  7^ 
Then  with  respect  to  the  set  y.  k  i  p+1,  9  Is  uniformly  robust  with 

order  N  and  there  Is  no  locking. 

Proof.  By  (4.7),  we  have  W(V*^)  =  Pp+j  j*  For  p  =  4,  this  gives  the 
Argyrls  triangle,  for  which  the  Interpolation  theory  (Theorem  6.1.1  of  [13]) 
shows  that 

C5.6)  Inf^  Ii*-X|L  -  S  ChPl|*|l  _ 

ph  2,H  p+2,n 

*  P+1.1 

Moreover,  the  proof  of  Theorem  6.1.1  of  [13]  may  easily  be  generalized  to 
show  that  (5.6)  holds  for  any  p  2  4.  The  result  follows  from  this.  (Note 
that  the  result  follows  Immediately  from  (5.5)  for  the  mesh  7^).  a 

Remark  5. 1.  In  [19],  a  more  restrictive  condition  on  the  mesh  Is  used  to 
prove  that  locking  does  not  occur  for  p  2  4.  However,  the  result  obtained 
there  Is  stronger  and  leads  to  divergence  stability  as  well,  as  discussed  in 
Remark  3. 1.  Ue  avoid  the  condition  on  the  mesh  precisely  because  our 
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Definition  3. 1  of  locking  is  different,  involving  only  the  displacements  and 
not  the  pressiires. 

So  far  we  have  discussed  the  triangular  mesh  7^  (and  7^).  Let  us  now 
discuss  the  mesh  7^  with  respect  to  the  extension  7  based  on  Q  (S).  We 
first  prove  some  auxiliary  lemmas. 


Lemma  5.3.  Let  =  {(x^.x^),  0  <  Xj.Xg  <  «>,  i  S  a  i  1,  be  covered  by 

the  mesh  Consider  the  set  G(x.,x.,)  *  G(x, )  =  x^^^  + 

2  p,— 1  2  12  11 

f(Xj)  where  f  is  a  polynomial  of  degree  S  p.  Then  for  1  »  1,2 


(5.7) 


Zi  =  inf^  _ 

n**€7j 


f  I 


IIG-Zll 


2,0“ 


1/2 


*  o(h**  . 


Proof.  Since  7^  is  a  uniform  mesh,  the  following  inverse  ineqviality  holds 
on  every  element  for  any  polynomial  v  e 


(5.8) 


l|v|ls,jjh  2S  Ch^'^llvllt^jj**  •  0  S  t  S  s. 


where  C  is  independent  of  v,  0  ,  t,  s,  h.  Suppose  for  i  ^  1,2 

Zj  =  o(hP"^) 

and  that  this  infimum  is  attained  by  Xj  €  .j-  Then  we  have,  using 

(5.8)  with  V  =  G-Xj  €  and  s  =  p+2-i,  t  =  2 


(I ' 

kj*‘€7j 


’*^”*i’*p+2-i,fl^ 


1/2 


S  =0(1) 


from  which 


(y  lx!-'  * 

'  *  sxr'-lo.n" 

2 


1/2 


=  0(1) 
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Implies 


0  over 


But  this  is  impossible,  since  Xj  «  ^^1-i  -i 


ax, 


p*2~l 


each  This  proves  (5.7). 

We  now  prove  the  following  lemma,  using  an  idea  from  the  proof  of 
Theorem  2  in  [7]. 


Lemma  5.2.  Let  «  {2^  be  as  in  Lemma  5.1  Euid  let  D(x^,X2) 


Then  for  i  =  1,2 


x^^^Xg.  Let 


W, 


(PL-..  . 


P+1.1 


inf  IID-zll-  _  *  o(hP‘^)  . 

X«Wi 


Proof.  We  use  the  standaird  notation 


(X  -  x)“  =  (X  -  x)‘ 
=  0 


Let  M  s  h 


-1 


and 


for  X  i  X 
for  X  <  X  . 

'  xJ)  . 


for  0  5  i,J  S  M-1,  let  (xj,  x^)  ■  (Ih,  Jh).  Then  any 


X  ®  Pp^j  can  be  represented  as  a  linear  combination  of  the  truncated 


powers 

(5.9) 


(x^-x');) 


0  5  i,J  S  M-1 


where  Pj»*Ij  ^  0  Pj  <Ij  ^  P  +  !•  Any  *  Qp  -j  ^  simileu'ly 

represented  by  taking  p^  S  p,  ‘Jj  ^  P*  I**  w®  Impose  the  constraints  p^  ^  2 
for  i  >  0  and  qj  k  2  for  J  >  0,  then  (5.9)  gives  us  precisely  all  z  € 
pjj  ,  and  qJJ  ,  respectively.  Consequently,  in  this  caise,  the  only  terms 
pj?  ,  of  the  form  (5.9)  with  p,  k  p  C8ui  be  (x -xj)^^^,  (x-xj)^  or 

P^  X  llT  1  + 


in 


p+i 

‘^p 


(x-x  )*^x  vdiile  those  satisfying  q,  2  p  must  be  of  the  form  (x-x^)^ 

1  2  J  2  ▼ 


(x-x;!)^  or  X  (x-x^)f. 
2  +  1  2  ♦ 
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For  H  >  0,  define  the  difference  operator  by 

2 

AjjgCXj.x^)  -  ^  [^jgCXj+la.  Xg+ma) 

l.n  «  0 

^diere  a  «  H/4.  Let  H  »  sh  vdiere  s  is  a  positive  integer.  Then  for  any 

X  «  p\,  ,  (of  the  fora  (5.9))  with  p,  ^  p  or  q,  2  p.  we  have 
p+1,1  1  j 

A„X  =  0. 

This  implies  that  AjjWj  c  j(^)  and  AjjWg  c  Q^_j  j(y2^’  Also,  we  see 
that 

AjjD  =  CxJ"^^  +  f(Xj) 


where  f(x^)  is  a  polynomial  of  degree  ^  p.  Hence,  if  we  fix  H  to  be 
sufficiently  small  and  let  s— »«  (i.e. ,  h— »0),  we  see  by  Lemma  5.1  that  for 
i  «  1,2 


inf  ||A„(D-x)||p  Q 
*«w,  ^  2,n, 


/2 


C  inf  ||x^ 

x€Qp+i-i ,1 


+1 


*  '■'*1'  -  * 


o(hP-‘) 


Noting  that 


(since  Just  gives  a  linear  combination)  proves  the  lemma.  □ 

Vfe  can  now  prove  theorems  for  rectangulau'  meshes  analogous  to  some 
theorems  proven  for  triangular  meshes. 


Theorem  5.5.  Let  the  extension  procedure  f  consist  of  the  h-version  on 
the  uniform  mesh  7^  for  elements  of  type  Q  with  p  >  1.  Then  with 

d  p 

1/2 

respect  to  the  solution  set  H.  ,  k  2  2,  7  shows  locking  of  order  N 

K,  V 

i.e.,  it  is  not  robust. 
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Proof.  The  proof  is  essentially  identical  to  that  of  Theorem  5.1.  We  now 

use  Theorem  1  (instead  of  Theorem  3)  of  [7]  to  characterize  the  density 

ofttespK.  :>  «,  uoj  By  this  theorem, 

2 

dense  in  C(R  )  only  when  s  2  4. 


is 

□ 


Theorem  5.6.  Let  the  extension  procedure  9  consist  of  the  h-version  on 
the  uniform  mesh  3^  for  elements  of  type  2p  with  p  i  2.  Then  with 
respect  to  the  solution  set  H.  ,  the  extension  9  is  uniformly  robust  with 

K|  V 

order  N  £qj,  k  ^  p  and  shows  locking  of  for  k  k  p+1. 


Proof.  By  the  results  in  [14],  the  optimal  rate  of  convergence  is  obtained 
when  tensor  product  splines  of  degree  k  2  in  each  variable  are  used. 
More  precisely,  we  have  for  ♦  e  k  k  p. 


,)=P"  .  .uQ;  2,Q  2,0 

p,0  Pfl  p,l 


Moreover,  by  Lemma  5.2,  this  is  the  best  rate  possible.  The  theorem 
follows,  using  Theorems  3.4  and  4.1.  □ 

We  now  consider  an  extension  procedure  based  on  7^  using  elements  of 
type  3p' . 

Theorem  5.7.  Let  the  extension  procedure  9  consist  of  the  h-version  on 

the  uniform  mesh  3,  for  elements  of  type  Q  '.  Then  with  respect  to  the 

c  P 

solution  set 

1  /P 

i)  For  p  »  1,  k  k  2,  f  shows  locking  of  0{N  )  and  hence  is 

non-robust . 

-1/2 

ii)  For  p  ■  2,  f  is  uniformly  robust  with  order  N  for  k  k  2 
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1/2 

and  shows  locking  of  order  N  for  k  k  3. 

Ill)  For  p  >  2,  f  Is  uniformly  robust  with  order 

for  k  >  p-t-2  and  shows  locking  of  order  N  for  k  k  p-t-1. 


Proof.  For  p  ^  2,  as  shown  In  Section  4, 


) 


). 


{fence,  by  Theorem  3.4,  the  results  for  locking  and  robustness  are  Identical 

for  ^  '  and  6  elements,  so  that  (1)  and  (11)  follow  from  Theorems 
P  P 

5.5  and  5.6,  respectively. 

By  (4.9)  and  Theorem  3.4,  we  see  that  for  p  >  2,  9  Is  robust  with 

uniform  order 


g(N).Sup,  ,  Inf 


By  Theorem  2  of  (7],  we  see  that 


Inf 


h 

p*i ,  1 


(Jj) 


''*‘^«c(n) 


5  (niP||*iij,(p*2)(jjj 


from  which,  using  the  Inverse  Inequality  and  the  fact  that  (q)  ^ 

^*****k+l,Q  obtain 

(5. 10)  i(N)  S  . 

Moreover,  by  Lemma  5.2,  this  Is  the  best  estimate  possible.  This  establishes 

(111).  □ 


Comparing  elements  of  types 
locking  but  that  the  locking  shown  by 


and 


Q  ' ,  we  see  that  both  show 
P 

elements  Is  twice  a  strong  as 


that  of  Q  . 

P 

Let  us  briefly  discuss  meshes  and  using  Q  elements.  For 

4  0  P 

the  mesh  5^,  the  statement  about  robustness  fo^N  exactly  the 
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sane  as  for  the  nesh  7^.  Nevertheless,  an  exact  assessment  of  the  locking 
is  not  presented  here. 

In  the  case  of  the  meshes  it  is  possible  to  generalize  the 

approach  of  the  construction  of  Bogner-Fox-Schmldt  rectangles  (see  [13])  to 
once  more  prove  robustness  of  0  (N  provided  p  k  4.  We  do  not 

discuss  the  case  p  <  4  here. 

Let  us  summarize  the  results  of  this  section: 


Type  of  mesh 

Type  of  element 

Order  of 
locking,  r 

f(N)=0(N^) 

Robustness 
order,  q 

g(N)=0(N*‘^) 

Uniform  triangular 

Pp.  P  »  1 

p  =  2 

p  =  3 

jBB 

mm 

Quasiuniform  triangular 

on 

r  =  0 

q  *  p/2 

Uniform  rectangular 

Ql.  Q^' 

r  =  1/2 

Qp.  P  i  2 

r  »  1/2 

q  »  (p-l)/2 

V'  P'2 

r  =  1/2 

q  =  1/2 

Qp'.  P  k  3 

r  =  1 

Quasiuniform  rectangular 

bb 

q  =  (p-l)/2 

Quasiuniform  parallelogram 

ma 

6.  The  D  and  h-o  versions. 

In  [23],  it  was  shown  that  the  p-verslon  (using  straight -sided 
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triangular  aeshes)  for  (3.1)  leads  to  a  robust  estimate  of  0[N  J,  c  > 

0  arbitrary,  vdien  the  solution  is  knovm  to  lie  in  i.e.,  the  order 

of  locking  is  at  most  N  .  The  loss  in  order  c  resulted  due  to  the 
definition  of  locking  employed  (see  Remark  3.1)  eind  is  directly  related  to 
the  locking  in  the  pressures.  Using  the  theorems  from  Section  3,  we  are  able 
to  obtain  an  optimal  robustness  estimate  and  prove  that  there  is  no  locking 
for  the  displacements  in  the  sense  of  Definition  3. 1.  This  holds  for 
parallelogram  elements  as  well.  Using  this  approach,  we  remove  the 
dependence  on  e  that  will  occur  in  the  robustness  estimate  for  the 
displeu:ements  if  the  results  from  [23]  are  used. 


Theorem  6. 1.  Let  the  extension  procedure  f  consist  of  the  p-version  using 

a  mesh  consisting  of  vriangles  or  parallelograms.  Then  with  solution  sets 

H.  ,  k  k  2,  f  is  free  of  locking  and  is  robust  with  uniform  order 
„-(k-l)/2 


Proof.  Using  the  results  of  [21]  and  [161,  we  have  for  the  p-version  with 
continuous  triangular  or  parallelogreun  straight -sided  elements  (h^ 
fixed),  for  ^  e  H‘'*^(n),  k  2  2, 

(6.1)  inf^ 

*  P.  1 

Since  c  W(^),  using  Theorem  3.4  gives  the  required  robustness  rate. 

The  absence  of  locking  follows  by  Theorem  4.2.  □ 

Using  the  results  of  [19],  it  can  be  shown  that  the  h-p  version  with 
quasiuniform  triaingular  meshes  results  in  locking  of  order  at  most  p^,  with 
the  rate  in  h  being  optimal,  provided  p  k  4.  We  are  able  to  prove  the 
following  theorem. 
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Theorea  6.2.  Let  the  extension  procedure  f  consist  of  the  h-p  version, 

using  quasiunifora  aeshes  consisting  of  triangles.  Let  the  solution  sets  be 

H.  with  k  ^  2.  Then  for  p  2  fflax(4.k-l},  9  is  free  of  locking  sind  is 
K,  V 

robust  with  uniform  order  N  more  precisely  h^ 

Proof.  Me  only  outline  the  proof  here.  Essentially,  the  idea  is  to  use 

(6.1)  together  with  a  standard  scaling  argiiment  to  show  that 

(6.2)  M 

provided  p  k  k-1.  For  details  see  [2],  where  this  scaling  airgument  h^LS  been 
used  to  prove  an  estimate  analogous  to  (6.2)  for  the  case  of  elements. 

Using  Theorems  3.4,  4.2  (or  Eqn.  (4.3))  completes  the  proof.  □ 


7.  Some  generalizations. 

So  far,  we  have  only  considered  the  two-dimensional  isotropic  case.  It 
turns  out  that  situations  analogous  to  Poisson  locking  ew'ise  in  more  general 
contexts  as  well.  For  example,  the  same  phenomenon  is  observed  in  the 
three-dimensional  Isotropic  case  when  v—>0.5.  Moreover,  one  can  consider 
simlleu'  situations  that  occur  in  the  anisotropic  case  (both  2-D  aind  3-D), 
when  locking  may  be  observed  due  to  the  introduction  of  a  constraint  on  the 
approximate  subspaces. 

Let  us  look  at  the  general  3-D  equations  of  anisotropic  eleisticity  given 


by 

(7. 1) 

(7.2) 

(7.3) 


A<r  >  g(u)  in  (2 

dlvo'  *  ?  in  n 

-♦  -»  _ 
<m  =  g  on  r 


where  <r  is  the  3x3  symmetric  stress  tensor,  and  c(u) 


the  strain  tensor. 
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A  is  a  fourth-order  tensor  known  as  the  compliance  tensor  which  depends  upon 
the  properties  of  the  material.  It  is  a  self-adjoint  lineau*  operator  acting 
on  the  six-dimensional  space  of  symmetric  3x3  tensors,  and  is  chau'acterized 
by  specifying  21  independent  elastic  moduli,  n  is  the  unit  normal  to  the 
boundary  8is  before  and  we  assume  (2.5)  again. 

The  standard  variational  form  of  (7. l)-(7.3)  is  obtained  under  the 
assumption  that  A  is  positive  definite,  in  which  case  (7.1)  may  be  solved 
for  <r  and  substituted  in  (7.2)  to  give  a  problem  involving  the  unknown  u 
alone.  For  many  important  materials,  however,  A  may  be  pxssitive 
semi-definite  and  singular  (or  close  to  singular).  In  this  case,  if  0  < 
i  Xg  . . .  ^  ^0  denote  its  eigenvalues  and  . .  .(r^  a  corresponding  set 

of  orthonormal  eigenvectors,  then  X^^  ®  0  (or  is  close  to  0).  The 
corresponding  constraint  on  u,  analogous  to  (1.1),  becomes 

(7. 4)  dlv  ((TjU)  =  0 

In  fact,  Poisson  locking  is  a  special  case  of  (7.4),  because  for  the  caise  of 
isotropic  materials  with  v  =  0.5,  it  can  be  shown  that  the  identity  matrix 
is  an  eigenvector  of  A  corresponding  to  the  eigenvalue  0. 

There  are  two  cases  which  must  be  distinguished  between.  (7.4)  defines 
a  singular  constraint  if  is  a  singular  tensor.  If  is  non-singular, 

then  (7.4)  is  called  a  non-singular  constraint.  Poisson  locking  is  obviously 
an  example  of  a  non-singular  constraint.  On  the  other  hsind,  a  material  that 
is  inextensible  in  the  direction  r  satisfies  (7.4)  with  a  singular  <r^ 
given  by 

(7.5)  (Tj  *  r^  r 

Both  singular  and  non-singular  constraints  will  lead  to  locking,  when 
the  material  is  such  that  A  is  (or  is  close  to  being)  singular.  Singular 
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constraints  like  (7.5)  are  analogous  to  the  locking  constraints  found  in 
shell  problems  and  are  generally  harder  to  develop  robust  methods  for.  See 
Section  3  of  [4],  where  we  have  discussed  the  robustness  of  various  methods 
for  the  problems  of  heat  transfer  through  highly  anisotropic  materials,  which 
has  the  same  character  as  locking  due  to  a  singular  constraint. 

Here,  we  are  interested  in  the  case  when  A  is  not  close  to  having  a 
singular  constraint  (this  is  characterized  precisely  in  [1])  but  is  close  to 
having  a  non-singular  constraint  like  (7.4).  Such  a  situation  has  been 
characterized  in  detail  for  the  case  of  3-D  orthotropic  materials  in  [1]. 
Denote  by  u.,  the  solution  to  (7.1)-(7.3)  when  the  compliance  tensor  A. 

Ai  Ai 

has  minimum  eigenvalue  >  0  and  let  the  limiting  constrained  case  be  A^. 
Then  by  Theorem  1.1  of  [1],  if  the  constraint  is  non-si ngulau*,  there  exists  a 
unique  limit  solution  u^  satisfying  (7.4)  with 


Lim 

*x.- 


"V“o'l.O  '  o 


It  is  expected  that  aui  analog  of  (2.12)  (with  A^  replacing  l-2i;)  will  be 
valid  in  this  case  as  well.  However,  such  a  result  is  not  currently 
available. 

Suppose  we  discretize  the  standard  vsiriational  formulation  of 
(7. l)-(7.3)  (obtained  by  first  solving  (7.1)  for  (r).  Then,  without 
establishing  (2.12),  Theorems  3.2  -  3.4  cannot  be  proved  as  stated.  However, 
it  is  still  instructive  to  look  at  the  limit  problem,  where  the  exact  emd 
approximate  solutions  must  satisfy  the  constraint  (7.4).  Suppose  we  are 
interested  in  exact  solutions  u  e  H,  where  H  cheuracterizes  the 
smoothness.  Let  Hq  be  the  subset  of  H  satisfying  (7.4)  and  Z(^)  the 
subset  of  satisfying  (7.4).  Then  the  limit  problem  hzis  the  optimal 

rate  of  convergence  Fq(N)  in  the  3^(0)  norm  if  and  only  if 
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(7.6) 


g(N)  -  §up  inf  IlS  -  w|L  -  S  CF-(N). 
o 

By  Theorems  2.2(A),  2.3(A)  of  [4],  (7.6)  Is  a  necessary  condition  for  5  to 
be  free  of  locking  with  respect  to  the  3^(0)  or  the  energy  norm. 

Vftien  0"^  is  non-singular,  we  can  derive  a  more  convenient  condition 
which  is  sufficient  for  (7.6)  to  hold  and  is  analogous  to  (3.11).  As  in  the 
two-dimensional  case,  we  have 

divz  =  0  ^  z  =  Curl 

Hence,  for  u  €  H^,  w  €  Z(V^),  we  have 

CjU  =  Curl  J,  o-jW  =  Curl  x 

where  I  €  W(H^),  w  €  W(V^).  Here, 

(7.7)  W(Y)  =  {J,  <r’Wl  J  6  Y)  . 


Then  we  have 


u  -  wllj  =  ll<r‘^  Curl  (?-Z)ili  jj  ^ 


so  that  the  following  is  sufficient  for  (7.6)  to  hold: 

(7.8)  Sup  inf  Il?-Z1L  o  i  CF  (N)  . 

;€W(H®)  z€W(v'‘)  ° 

As  em  application  of  (7.8),  let  us  look  at  the  case  that  {^}  consists 
of  continuous  piecewise  polynomials  of  fixed  degree  p  on  a  quasiuniform 
sequence  of  meshes  consisting  of  tetrahedra.  Then  for  any  non-singular  o-^ . 
will  be  bijective  on  so  that  from  (7.7),  we  see  that 

w(\^)  =  v(^) » pj; 

P  P+1.1 

analogously  to  (4.7).  By  [24],  the  minimum  q  required  to  obtain  an  optimal 
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rate  of  convergence  for  the  h-verslon  using  pIJ  ,  Is  q  »  9,  so  that  (7.6) 

* 

will  be  satisfied  whenever  with  p  ^  8  and  the  h-verslon  for  the 

limit  problem  will  show  optimal  convergence.  We  expect  once  again  (similar 
to  the  results  in  Section  5)  that  the  use  of  lower  p  or  the  use  of 
peirallelepiped  elements  will  result  in  locking. 
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